General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




X-692-70-4S2 

■ SPRINT 


im Iti to $*/ $ 6 

ON THE NATURE AND ORIGIN OF 
DIRECTIONAL DISCONTINUITIES 


L. F. BURLAGA 


DECEMBER 1970 




k 9 Si A 

>> cm 

k ^ 




- GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 

N 71- 167 5 7 

(ACCESSION. NUMBER) ’ 

/y _- — .< 


(THRU! 

r\ o 


(PAGES) 


1 M- (,<4^0 


=f ~t~ — ' - , 

(NASA CR OR TMX OR AD NUMBER) 


(CODE) 

(CATEGORY) 


P . 


ON THE NATURE AND ORIGIN OF 


DIRECTIONAL DISCONT"> i* , '"IES 


Leonard F. Burlaga 

Laboratory for Extraterrestrial Physics 
NASA -Goddard Space Flight Center 
Greenbelt, Maryland 20771 


November 1970 


Abstract 


Less than 257„ of the "directional discontinuities" (directional 
changes >30° in less than 30 sec) identified by Burlaga in the 
interplanetary magnetic field data of Ness from Pioneer 6 are Alfven shocks 
( = rotational discontinuities). Thus, most of these directional 
discontinuities are tangential. The Pioneer 6 magnetic field and plasma 
data are not consistent with the hypothesis that directional discontinuities 
originate primarily by the collision of fast streams with slower plasma, 
since the rate of occurance of discontinuities in regions of increasing 
bulk speed is only a little higher than the rate elsewhere. The rate at 
which discontinuities are seen at .82 AU, is somewhat lower than at 1 AU, 
but the difference may be due to the lower data quality and increase in the 
number of data gaps when the spacecraft is far from the earth. The 
distributions of the change in direction of the field across the 
discontinuity and of the "normals" (B^xB^/B^xB^) are essentially the same* 
at «.8, «.9 and «1.0 AU. This suggests that most discontinuities originate 
within .8 AU and do not evolve appreciably between .8 and 1.0 AU. 
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INTRODUCTION 

Discontinuities in the interplanetary magnetic field on a scale of .01 
AU were discovered by Ness et al. ( 1966 ). Statistical studies by Siscoe 
et al . ( 1968 ) and Burlaga ( 1969 a) showed that these discontinuities pass 
a spacecraft at the rate of«l/hr. Since fast or slow shocks occur much 
less frequently, and since contact discontinuities show no change in B, the 

discontinuities found by Ness must be predominantly tangential discontinuities, 
Alfven shocks, or a mixture of both (see Colburn and Sonett ( 1966 ) for a 
discussion of the types of hydromagnetic discontinuities Burlaga ( 1968 ) an d 
Burlaga and Ness (1969) used magnetic field and plasma data to show that 

some of the discontinuities are tangential. Siscoe et al. (1968) suggested 
that, in fact, essentially all interplanetary discontinuities are tangential 
discontinuities. On the other hand, Smith et al. (1970), Belcher and Davis 

(private communication ), and Turner and Siscoe (1970) suggested that most 
'discontinuities ' 1 are Alfven shocks. Thus, the relative number of tangential 
discontinuities and Alfven shocks is controversial. 

In view of the great complexity of interplanetary magnetic field 
changes, and noting the possibility that the discontinuities found by 
Ness might be a mixture of 2 types, Burlaga ( 1969 a) introduced the 

term "directional discontinuity" to denote a discontinuity in which the 

o 

magnetic field direction, uj, changes by >30 in < 30 sec. Most of the 
discontinuities identified by Ness al. (1966) and Siscoe et al. (1968) are 
directional discontinuities. It is not clear whether or not the discontinuities 
referred to be Smith et al., Belcher and Davis, and Turner and Siscoe are 
directional discontinuities. 
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This paper has two aims: l) to determine the ratio of Alfven 

shocks to tangential discontinuities in the set of directional discontinuities 
obtained by Burlaga (1969a) from Ness’s Pioneer 6 magnetic field data for 
the period Dec. 16, 1965 to Jan. 5, 1966; and 2) to discuss some hypotheses 
concerning the origin of these discontinuities, a problem that has not 
been studied experimentally until now. 
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II. Nature of Directional Discontinuities 

There are a few ways to distinguish between a tangential 
discontinuity and an Alfven shock, but only one could be applied to 
the data under consideration. If a discontinuity is an Alfven shock, 
then 



(1) 


here 

a = (i , (2) 

is a measure of the anisotropy of the velocity distribution, 

\/ is the bulk velocity, B is the magnetic field, p is the plasma density, 
and P |( and P^« are the total thermal pressure (ions plus electrons) parallel 
to or perpendicular to EL In general, |b^| ^ |fipj 
in which case P ^ ^ ^ (Hudson, 1970). 

If a discontinuity is tangential, then AV=Vp - can have any value 
and is independent of B (but very large AV may give rise to the Kelvin- 
Helmholtz instability; Parker (1963)» Burlaga (1969b)). Thus, one way to 
set an upper limit on the number of Alfven shocks is to determine the 
number of discontinuities which satisfy (l). This givs an upper limit because 
some tangential discontinuities might also satisfy (l) by coincidence. 

There is one complication with the metliod just described: The 

anistropy was not measured for the discontinuities under consideration, so 


A cannot be determined for each discontinuity. However, the value can 
be reliably estimaled as follows. Rrom Hundhausen et. al. (19?0), 
(T maY /T T T n -y 1 ) median = 1.9» since the thermal anisotropy is 

field aligned, this is T /f l for the protons. Measurements for the 


■ 


5 


electrons are not cl ear cut, but Montgomery et al. (1968) find that 
nearly always (T e /T e . ) <1.2. Thus , we shall take 

T e ||/T e i^l.l. Hundhausen et al. (1970) give 6.9xlc/ l °K and 6.6 cm ^ for 
the median proton temperature and density, respectively, and the electron 
temperature is«1.5xl0^°{< (see Montgomery et. al. (1968), and the brief 
summary in Burlaga and Ogilvie (1970 )). The magnetic field intensity is 
typically n ’6y. With these numbers, one finds that A =(.9 ± .1). 

We shall consider a subset of the discontinuities listed in Burlaga 
(1969a), those for which accurate plasma data are available from the MIT 
plasma instrument on Pioneer 6. The subset consists of 200 directional 
discontinuities in the period Dec. 17, 1965 to Jan. k, 19 66. Bi and Bo 
were obtained as described in Burlaga (1969a). The parameters n^, V 
n£ and were obtained 

from the spectra just before and after the discontinuity, 
except when the discontinuity might have occured while a spectrum was being 
recorded, in which case the preceding (or following) spectrum was used. A 
discussion of the plasma instrument may be found in Lazarus et al. (1966). 


The vectors B and V were resolved in BTN components, where R is radial 

/> 

away from the sun, T is in the ecliptic plane and points in the direction 
of the earth’s motion, and n is perpendicular to the ecliptic. The 
most accurately measured velocity component is v^, the least accurate is 


N* 


In these coordinates, (l) may be written 


<V V 


± A > 


i - R,T,N, 




(3) 


where 0 = 21.8 (-li - — £i) (Hi) ^ (1+) 

1 nj n 2 lt n 

and n is the number density. 

All of the quantities on the LHS of (3) are known. The question at 


6 


issue may now be stated as follows: "What fraction of the directional 
discontinuities satisfy ( 3 ) with A w *9 + .1?" 

Figure 1 shows the distribution of AV./q. for the 200 directional 

i 

discontinuities. If the discontinuities were mostly Alfven shocks the 
distributions would be peaked at +.9 of -,9.The observations in Figure 1 
are not distributed in this way. Rather , the distributions are all 
peaked at zero, as one might expect for tangential discontinuities. 

Clearly, then, most of the directional discontinuities are not Alfven 
shocks, and are thus probably tangential discontinuities. 

A quantitative upper limit on the fraction of Alfven snocks may be 
obtained as follows. Figure 2 shows AV versus 0^. If the discontinuities 
were all Alfven shocks, the points would scatter about the line AV=,9 . 
Suppose that all of the points above the line|Av( = ,9Q in Figure 2 
are rotational discontinuities and that an equal number (26) of the points 
below the line are rotational discontinuities. Since there are 200 
directional discontinuities in the sample, this implies that <2%% of the 
directional discontinuities are rotational .. This is clearly an upper 
limit, since some of the points above the line in the region Q<10 are 
are almost certainly tangential discontinuities. 

Our sample of discontinuities includes both a very quiet time 
(v <350 km/sec) and a very disturbed period (V<600 km/sec.) The 
distribution of the discontinuities on the most magnetically disturbed 
day in our sample £hown by the crosses in Figure 2) is not significantly 
different from the general distribution. 
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III. Origin of Directional Discontinuities . 

We shall examine two hypotheses: 

1) . Discontinuities are caused by the impact of fast streams with 
slower plasma near 1 AU. 

2) . Discontinuities evolve by some process operating between 
,8 AU and 1 AU. 

Impact Hypothesis . Figure 3 shows directional discontinuities for 
the period December 16 , 1965, to January 5* 1966, as given by Burlaga 
(1969a) , together with an indication of the times when the bulk speed 
was increasing for an hour or more (shaded areas). Cross hatched areas 
show times for which high bit-rate magnetic field data are not available, 
and the double-arrowed lines above the panels indicate times for which 
the plasma data were not available. Clearly, there is no strong tendency 
for directional discontinuities to cluster in regions where the bulk 
speed increases r so the faster streams cannot be the principal cause of 
directional discontinuities . However, the rate of occurance of discontin** 
uitiee in the gradient regions, (l.5±. 2 )/hr, is possibly somewhat higher 
than outside these regions, (l.04>l)/ir. 

Evolution Hypothesis . This implies that the statistical properties 
should change with distance between .8 AU and 1 AU. We shall consider 
3 distances, .82 AU (perigee), .91 AU, (where the trajectory crossed the 
earth-sun line) and .98 AU. The corresponding time intervals are 2000 UT 
April 29 through 1900 UT June 5, 1966, Feb. 25 through March 8, 1966, and 
Dec. 16, 1965 through Jan. 5, 1966, respectively. 

The statistical properties that we shall consider axe a) the rate 
at which discontinuities pass the spacecraft (this is a measure of their 
"density” , since the solar wind speed does not change appreciably between 




8 


.8 and 1 AU) , b) the distribution of u), tho change in the magnetic field 
direction across a discontinuity, and c) tie distribution of "normals", 
ft ■ B^xB 0 . The statistics are lees accurate at .82 and .91 AU because of 
data gaps, lower telemetry rates and noise. 

Subjectively selecting hours with "good" data for the periods of 

interest gives .7 discontinuities/lir at .82 AU, ,8 at .91 AU, and 1.1 at 

,98 AU. The higher frequency of discontinuities near the earth is at least 
partly, and possibly entirely, due to the better and more continuous data 
near the earth. As a measure of the activity during these periods, we 
note that was 12.0, 8.6 and 12.3 for Dec. 16-Jan. 3, Feb. 25-Mar. 8, 
and April 29- June 5, respectively. In view of the (weak) correlation 
between and the density of discontinuities (Burlaga, 1969a), one might 
expect the first and last periods to be comparable, and the frequency of 
discontinuities for the intermediate period to be relatively low. We 
conclude that the density of directional discontinuities at .82 AU is 

possibly 35% less than at .98 AU, but the difference may be due mostly to 

the higher quality data near 1 AU. 

The distributions ofou are shown in Figure 1+. There is no significant 
difference between the 3 distributions, indicating little change in 
discontinuities as they move from .82 AU to 1 AU. 

The distributions of "normals" (B xB o ) are shown in Figure 5 . These 
distributions were computed by putting all of the vectors in their proper 
orientations at the origin of a urit sphere and computing the density of 
the points which were determined by the intersections of these -vectors with 
the unit spheres. The left side of Figure 5 shows the density integrated 
over as a function of 9. (9 is the angle between the "normal" and its 
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projection in the ecliptic plane, and 6 is the angle between the projection 
and the earth-sun line). The right side of Figure 5 shows the density, 
integrated over 0, as a function of 6 . The radius of each sector is 
proportional to the number of normals per unit area on a unit sphere, 
for that interval. There is no significant difference between the 0 
distributions for .82 AU, .91 AU, and .98 AU. The 6 distributions at .82 
and .91 AU are essentially the same, but appear to differ from that at 
.98 AU. The difference is probably not statistically significant, however, 
since there are only 25 to 30 discontinuities in the largest sectors. 

We conclude that most if not all directional discontinuities originate 
within .82 AU and do not change appreciably between .82 and .98 AU. 
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Figure 1 


Figure 2 


Figure 3 


Figure 1+ 


FIGURE CAPTIONS 

Distribution of AV /Q,^. If the discontinuities were Alfven 
shocks the distribution would be peaked near +.9. 

The change in bulk speed across directional discontinuities 
versus the corresponding change in Qf=21. 8 jB^/nj-B^/n^l (^I)^ 2 
where B R is the radial component of B. (See text) If all 
directional discontinuities were Alfven shock, the points would 
scatter about the lines at +.9 in this figure. 

Relation between directional discontinuities and increasing 
bulk speed. The discontinuities are shown by the vertical 
lines whose position gives the time they were observed 
at Pioneer 6 and whose height shows the angle through 
which B changes across the discontinuity. The shaded 
regions indicate times when the bulk speed was increasing. 

There is no clear association between the shaded regions 
and discontinuities, indicating that fast streams are 
not the principal cause of directional discontinuities. 

The distribution of u) for discontinuities observed at 
.98 AU, .91 AU, and .82 AU. No significant difference 
between the distributions can be seen, supporting the 
b othesis that directional discontinuities do not evolve 
between .82 and .98 AU. 

Distribution of the direction of the vectors B_xB_ for the 

—2 

directional discontinuities at .82, .91, and .98 AU. 


Figure 5 
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Figure 4 




